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In this paper we compare systematically the two-component Dual Parton Model (DPM) event generator Dpm- 
jet-III to d-Au and p-p data from RHIC. In this process we are able to improve the model. The need for fusion 
of chains and a recalibration of the model to obtain collision scaling in h-A and d-A collisions was found already 
in previous comparisons. Here, comparing to transverse momentum distributions of identified charged hadrons 
we find also the need to modify the transverse momentum distributions in the decay of hadronic strings, the basic 
building blocks of the model on soft hadronic collisions. 



1. Introduction 

Hadronic collisions at high energies involve the 
production of particles with low transverse mo- 
menta, the so-called soft multiparticle produc- 
tion. The theoretical tools available at present 
are not sufficient to understand this feature 
from QCD alone and phenomenological models 
are typically applied in addition to perturbative 
QCD. The Dual Parton Model (DPM) [T] is such 
a phenomenological model and its fundamental 
ideas are presently the basis of many of the Monte 
Carlo implementations of soft interactions. 

The soft component of the Dual Parton Model 
contains some features to be determined from 
multi-particle production data. This is so for the 
fragmentation of chains. But it also concerns the 
general event structure. Already in Ref. we 
have used data from RHIC to tune some of the 
properties of the soft particle production model 
contained in Dpmjet-III. This process will be 
continued here using mainly RHIC date on d-Au 
and p -p collisions. 

The Dpmjet-III code system |3I4| is a Monte 
Carlo event generator implementing Gribov- 

* Based on a contribution to a workshop "QCD at Cosmic 
Energies, Erice, Italie, Aug.29 to Sept. 5 2004 



Glauber theory for collisions involving nuclei, for 
all elementary collisions it uses the DPM and 
(LO) perturbative QCD as implemented in Pho- 
jetI.12 |SEj. 

Dpmjet-III is unique in its wide range of 
application simulating hadron-hadron, hadron- 
nucleus, nucleus- nucleus, photon-hadron, photon- 
photon and photon-nucleus interactions from a 
few GeV up to cosmic ray energies. 

The properties of Dpmjet-III were recently 
summarized in Ref. [2]. For a more dataillcd 
description of Dpmjet-III we refer to this paper 
and the literature quoted there. 

2. Comparing the original Dpmjet III 
with RHIC data 

In Ref. [2] we found an excellent agreement of 
Dpmjet-III up to transverse momenta of about 
p± = 10 GeV/c to 7T° transverse momentum dis- 
tribution in p-p collisions at yfs — 200 GeV from 
PHENIX 0. 

In the same paper, comparing Dpmjet-III 
with multiplicity data from RHIC in central Au- 
Au collisions we found a problem which let us 
introduce percolation and chain fusion into Dp- 
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mjet-III. The groups at Lisboa [819110) and San- 
tiago de Compostela ^5 were the first to point 
out, that the multiplicities measured at RHIC are 
significantly lower than predicted by conventional 
multi-string models. A new process is needed to 
lower the multiplicity in situations with a very 
high density of produced hadrons like in central 
nucleus-nucleus collisions. The percolation pro- 
cess, which leads with increasing density to more 
and more fusion of strings, is one such mechanism 

nam 

3. Pseudorapidity distributions in d-Au 
and p — p collisions and Dpmjet-III 

Pseudorapidity distribution of charged hadrons 
produced in minimum bias y/s = 200 GeV 
d-Au collisions were measured at RHIC by 
the PHOBOS-CollaborationO and by the 
BRAHMS CollaborationQ2| . In Fig [J we com- 
pare the preliminary PHOBOS data and the pre- 
liminary BRAHMS data to Dpmjet-III calcula- 
tions. Using Dpmjet-III without fusion of chains 
we find the Dpmjet distribution above the exper- 
imental data outside the systematic errors. Using 
Dpmjet-III with fusion of chains we find the Dp- 
mjet distribution within the systematic errors of 
both experiments. 

In [TB| the PHOBOS Collaboration also 
presents preliminary data on the pseudorapidity 
distribution of charged hadrons in p-p collisions 
at *Js — 200 GeV. In FigH] we compare Dpm- 
jet-III also with these data and find an excellent 
agreement. 

Pseudorapidity distributions of charged ha- 
drons in d-Au collisions with specified centralities 
were given by the BRAHMS-Collaboration[T5] 
and by the PHOBOS-Collaboration 17 . In Fig[2 
we compare the preliminary PHOBOS-data with 
results from Dpmjet. In the most central colli- 
sions we find an excellent agreement. In Dpmjet- 
III we know the impact parameter of the collision 
and we determine the centrality according to the 
impact parameter distribution. This might give 
results slightly different from the PHOBOS cen- 
trality determination, also the the number of par- 
ticipants determined by PHOBOS and obtained 
from Dpmjet-III differs for centralities between 
20 and 80 % , see Table 1. 



DPMJET (ns) d-Au 200 GeV m.b. 
PHOBOS d-Au 200 GeV m.b. 
BRAHMS d-Au 200 GeV m.b. 
DPMJET cha p-p 200 GeV m.b. 
PHOBOS cha p-p 200 GeV m. b. 
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Figure 1. Pseudorapidity distribution of charged 
hadrons produced in minimum bias i/s = 200 
GeV d-Au and p-p collisions. The results 
of Dpmjet with fusion of chains are com- 
pared to preliminary experimental data from the 
BRAHMS-CollaborationQ2] and the PHOBOS- 
Collaboration [l"lil6) . At some pseudorapidity 
values the systematic PHOBOS-errors as es- 
timated by the experimental collaboration are 
given. 



Table 1. Numbers of participants in d-Au colli- 
sions with different centralities. PHOBOS num- 
bers taken from ^7] 



Centrality % 


source 


N part 


iVpart(Au) 


AWt(d) 


0-20 % 


PHOBOS 


15.5 


13.5 


2.0 




Dpmjet 


15.4 


13.5 


2.0 


20-40 % 


PHOBOS 


10.8 


8.9 


1.9 




Dpmjet 


11.7 


9.8 


1.94 


40-60 % 


PHOBOS 


7.2 


5.4 


1.7 




Dpmjet 


7.8 


6.0 


1.76 


60-80 % 


PHOBOS 


4.2 


2.9 


1.4 




Dpmjet 


4.7 


3.24 


1.41 


80-100 % 


PHOBOS 


2.7 


1.6 


1.1 




Dpmjet 


2.9 


1.8 


1.1 
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Figure 2. Pseudorapidity distribution of charged 
hadrons produced in non-single-diffractive (ns) 
y/s = 200 GeV d-Au collisions with different cen- 
tralizes. The results of Dpmjet with fusion of 
chains are compared to preliminary experimental 
data from the PHOBOS-Collaboration[r7| . 

4. Transverse momentum and transverse 
energy distributions of charged hadrons 
in d-Au collisions and Dpmjet-III 

Let us repeat an important Dpmjet modifica- 
tion introduced already in our previous paper [5] . 
Several RHIC experiments (see for instance |18|1 
found in d-Au collisions at large p± a nearly per- 
fect collision scaling for ir° production. Collision 
scaling means Raa ~ 1.0 where the Raa ratios 
are defined as follows: 



Raa — 



dp±dr/ 



N 



A-A 



N, 



d 2 



A-A 

binary dp±dr} 



N N- 



■N 



(1) 



Here N^~^, y is the number of binary Glauber col- 
lisions in the nucleus-nucleus collision A-A. 

Dpmjet-III in its original form gave for ir° pro- 
duction in d+Au collisions strong deviations from 
collision scaling (Raa ~ 0-5 at large p±). The 
reason is understood^- Using a modified itera- 
tion procedure it was possible in |2] to obtain a 
nearly perfect collision scaling. 

Having shown the collision scaling of the 
model in |2| we will here compare the model 
directly with the measured transverse mo- 
mentum and transverse energy distributions. 
Transverse momentum distributions in mini- 



mum bias d-Au collisions were measured by the 
PHENIX-Collaboration[]J| and by the STAR- 
Collaboration 20 . In Fig|3] we compare the pre- 
liminary PHENIX-data for positively charged 
hadrons with Dpmjet-III as described so far. 
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Figure 3. Transverse momentum distributions of 
positively charged hadrons in minimum bias d- 
Au collisions. Compared are the preliminary data 
from the PHENIX-Collaboration[TU| to the re- 
sults of the original Dpmjet-III. 

For charged Pion production we find in Fig|3] 
quite a good agreement between Dpmjet-III and 
the preliminary PHENIX-data. This agreement 
is similar to the good agreement found already 
in Ref. [2] for ir° production in p-p and d-Au 
collisions. However, it seems, that the Dpmjet- 
III model fails to describe the p± distributions 
of heavier hadrons like Kaons and protons. The 
prediction is slightly too steep. This is the first 
time that we are able to compare the Dpmjet- 
III model to p± distributions of identified hadrons 
like K + , K _ , protons and antiprotons. In the past 
usually only data for all charged hadrons were 
available. Therefore, it is not surprising to find 
problems. 

The reason for this disagreement is the follow- 
ing: 

Dpmjet uses the Pythia code fo r the 
fragmentation of all chains (hard, QCD based 



chains as well as soft chains). Pythia selects 
in the PYPTDI function the transverse momenta 
of (for instance) q — q pairs (to become hadrons 
finally) from independent Gaussians in the p± x 
and p± y components of p±. Pythia was exten- 
sively tested in processes with a significant hard 
component, for the fragmentation of hard chains 
this choice is reasonable, but the fragmentation 
of the soft component is not well determined. 
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Figure 4. Transverse momentum distributions 
of positively charged hadrons in minimum bias 
d-Au collisions. Compared are the preliminary 
data from the PHENIX-Collaboration[T3| to the 
results of Dpmjet-III with modified transverse 
momentum distribution in hadronic soft chain de- 
cay as described in the paper. 

Otherwise, it is known since very many years 
(given already in 1968 in Ref. [22], but this might 
even not be the first Reference) that the trans- 
verse momentum distribution for hadrons with 
mass m in soft hadronic collisions is well described 
by an exponential in the transverse energy 



exp(-TOi/T), 



(2) 



We here chose a similar parametrization and re- 
place the Gaussian in PYPTDI of Pythia by the 
distribution 



After such a change all Pythia parameter which 
relate to the fragmentation have to be reopti- 
mized. 

With this modification we compare Dpmjet- 
III again to the PHENIX-transverse momentum 
distributions m FigQ] for positively charged 
hadrons and in Fig[S] for negatively charged 
hadrons. To be definitive: we use Dpmjet-III 
with (i) chain fusion (ii) anomalous baryon 
stopping as described in Ref. [23], (hi) the mod- 
ified iteration procedure to obtain collision scal- 
ing [2] and (iv) the changed transverse momen- 
tum distribution in soft chain decay as described 
above. We find in both Figures a much improved 
agreement to the preliminary PHENIX-data . 
At sufficiently large p± the p± distributions of Pi- 
ons, Kaons and baryons have the same exponen- 
tial slope. 
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exp(-Jpl + ml/ a), m x = OMGeV/c. (3) 
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Figure 5. Transverse momentum distributions 
of negatively charged hadrons in minimum bias 
d-Au collisions. Compared are the preliminary 
data from the PrIENIX-Collaboration[j"5] to the 
results of Dpmjet-III with modified transverse 
momentum distribution in hadronic soft chain de- 
cay as described in the paper. 

As a further check we present in Figs. and 
the comparison of the modified Dpmjet-III 
with tansverse energy distributions of identified 
charged hadrons as measured by the PH0B0S- 
Collaboration 24 . Before the modification there 
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were discrepancies between the Kaon and baryon 
distributions similar to the ones in Fig|3l The 
agreement in Figs. H3 and is now much better. 
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Figure 6. Transverse energy distributions of pos- 
itively charged hadrons in minimum bias d-Au 
collisions. Compared are the preliminary data 
from the PHOBOS-Collaborationgl to the re- 
sults of Dpmjet-III with modified transverse 
momentum distribution in hadronic soft chain de- 
cay as described in the paper. 

Finally in FiglHlwe compare the modified Dp- 
mjet-III with the transverse momentum distri- 
bution of all charged hadrons as measured by the 
PHENIX-CollaborationP) . We find a reason- 
able agreement. 

5. Transverse momentum distributions of 
charged hadrons at different pseudora- 
pidities 

Transverse momentum distributions of all 
charged hadrons and of all negatively charged 
hadrons at different pseudorapidities were mea- 
sured by the BRAHMS-CollaborationgSj in d- 
Au and p-p collisions. The interesting point in 
these measurements is the gradual change of the 
shape of the p± distribution with changing pseu- 
dorapidity and the corresponding evolution of the 
nuclear modification factors. In Figs. IH1 and 1X771 we 
compare the modified Dpmjet-III as described 
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Figure 7. Transverse energy distributions of neg- 
atively charged hadrons in minimum bias d-Au 
collisions. Compared are the preliminary data 
from the PHOBOS-Collaboration[23 to the re- 
sults of Dpmjet-III with modified transverse 
momentum distribution in hadronic soft chain de- 
cay as described in the paper. 

in this paper with the preliminary d-Au and p-p 
data of the BRAHMS-Collaboration[2Hj at pseu- 
dorapidities of 0, 1, 2.2 and 3.2. In the Dpmjet- 
III colculations we are not able to simulate the 
quite complicated pseudorapidity acceptances as 
given in [2H], instead we use in all cases a pseudo- 
rapidity band of width 0.2 centered at the nom- 
inal pseudorapidity. This difference might be re- 
sponsable for part of the differences we find be- 
tween the model and the data. There is no essen- 
tial difference between the agreements for d-Au 
and p-p. The model follows the change of the 
shapes of the data with changing pseudorapidity, 
but the agreement is not perfect. 

6. p — p collisions in Dpmjet-III 

We are here mainly concerned with d-Au colli- 
sions, but the changes of Dpmjet described in the 
last Section will also change hadron production 
in p-p collisions. Therefore, let us shortly discuss 
p-p collisions at the energy of the RHIC collider. 

In FigQ] we did already compare Dpmjet- 
III with the preliminary data of the PHOBOS 
Collaboration^) for the pseudorapidity distribu- 
tion of charged hadrons in p-p collisions at y/s = 
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Figure 8. Transverse momentum distributions of 
all charged hadrons in minimum bias d-Au colli- 
sions. Compared are the preliminary data from 
the PHENIX-CollaborationfTB) to the results of 
Dpmjet-III with modified transverse momentum 
distribution in hadronic soft chain decay as de- 
scribed in the paper. 



200 GeV. We found an excellent agreement simi- 
lar to the agreement found in earlier comparisons 
with data from the CERN-SPS and the TEVA- 
TRON colliders. 

In Figs. [5] and El we compared the modi- 
fied Dpmjet-III as described in this paper also 
with the p-p preliminary data of the BRAHMS- 
Collaboration 25 on charged hadron transverse 
momentum distributions at pseudorapidities of 0, 
1, 2.2 and 3.2. 

Using Dpmjet-III with modified transverse 
momentum distributions in hadronic soft chain 
decay in Figs El an d El we find like in d-Au 
collisions an improved agreement between Dp- 
mjet-III and the preliminary data on tt + , ir~ 
K + , K _ , proton and antiproton transverse mo- 
mentum distributions in p-p collisions from the 

STAR-Collaboration[2Dl ■ 

Transverse momenum distributions of neutral 
strange hadrons in minimum bias p-p collisions 
were measured by the STAR-Collaboration|26j . 
In FigEHwe compare Dpmjet-III with modified 
transverse momentum distributions in hadronic 
soft chain decay to these preliminary STAR data. 
The ageement is satisfactory. 



7. Summary 

The data obtained at RHIC are extremely use- 
ful to improve hadron production models like 
Dpmjet-III. 

Of particular importance are data on hadron 
production in p-p collisions, d-Au collisions and 
peripheral Au-Au collisions, in all of these colli- 
sions (unlike central Au-Au collisions) we do not 
expect any change in the reaction mechanism, 
which might not be accommodated into the mech- 
anisms as implemented in Dpmjet-III. 

Indeed, comparing Dpmjet-III to RHIC data 
we find three important corrections to be ap- 
plied to Dpmjet-III, which otherwise do not 
completely change the independent chain frag- 
mentation model: (i) Percolation and fusion of 
chains, the data from RHIC allow to determine 
the amount of percolation to be implemented 
into Dpmjet-III. (ii) Collision scaling of large p± 
hadron production in d-Au collisions: The data 
indicate that we have to change the iteration pro- 
cedure (of the selection of all soft and hard chains 
in nuclear collisions) in such a way, that colli- 
sion scaling is obtained, (iii) Replacing in soft 
hadronic collisions the Gaussian transverse mo- 
mentum distribution contained in the Jetset- 
Pythia code by an exponential distribution. 
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Figure 9. Transverse momentum distributions of 
all charged hadrons at r] crn = 0, at r] cm = 1, and 
at rjcm — 2.2 in minimum bias p-p and d-Au 
collisions. See caption FigfTUl 
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Figure 10. Transverse momentum distributions 
of all negatively charged hadrons at r\ crn = 
3.2 in minimum bias p-p and d-Au collisions. 
Compared are the preliminary data from the 
BRAHMS-Collaboration[55] to the results of Dp- 
mjet-III with modified transverse momentum 
distribution in hadronic soft chain decay as de- 
scribed in the paper. 
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Figure 11. Transverse momenum distributions of 
positively charged hadrons in minimum bias p- 
p collisions. Compared are the preliminary data 
from the STAR-Collaborationj20| to the results 
of the Dpmjet-III with modified transverse mo- 
mentum distributions in hadronic soft chain de- 
cay as decribed in this paper. 
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Figure 12. Transverse momenum distributions of 
negatively charged hadrons in minimum bias p- 
p collisions. Compared are the preliminary data 
from the STAR-Collaboration|2T)| to the results 
of the Dpmjet-III with modified transverse mo- 
mentum distributions in hadronic soft chain de- 
cay as decribed in this paper. 
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Figure 13. Transverse momenum distributions of 
neutral strange hadrons in minimum bias p-p col- 
lisions. Compared are the preliminary data from 
the STAR-Collaboration[21! to the results of the 
Dpmjet-III with modified transverse momentum 
distributions in hadronic soft chain decay as de- 
cribed in this paper. 



